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Abstract. We construct Green's function for the quantum degenerate parametric oscillator in 
terms of standard solutions of Ince's equation in a framework of a general approach to harmonic 
' oscillators. Exact time-dependent wave functions and their connections with dynamical invariants 

^SJ I and SU (1,1) group are also discussed. 
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1. The Degenerate Parametric Amplifier Hamiltonian 



p We use standard annihilation and creation operators in the coordinate representation [H], |20 
given by 



with p = —ihd/dx and oa^ — a^a = 1 throughout this Letter. The following Hamiltonian 

*> I H{t) = ^{aa^ + a^a)-^ (e^-* a' + e'^-* (at)') (1.2) 

^ ■ of the time-dependent Schrodinger equation 

tn^=H(t)ij (1.3) 

^ describes the process of degenerate parametric amplification in quantum optics. The first term 
corresponds to the self-energy of the oscillator representing the mode of interest, the second term 
describes the coupling of the classical pump to that mode (with the phenomenological constant A), 
giving rise to the parametric amplification process (the phase of the pump is taken to be zero at 
i = for convenience). Takahasi [15] wrote a detailed paper on this subject as far as in 1965 (see 
.=1 ; also [21], [33], [35], [36], [37], [38], [39], [H], W\ and references therein). The same Hamiltonian 
had also been considered by Angelow and Trifonov in the mid 1990s in order to describe the light 
propagation in the nonlinear Ti : LiNbO^ anisotropic waveguide (details can be found in [1] and 
[2]). Another form of this variable quadratic Hamiltonian is given by 

H = -L(^i + Acos(2u;t)^p2^^(^l-^cos(2u;t)^a;2 (1.4) 

A h d 

+ -sm(2ut) (px + xp] , P=-7^- 
2 t ox 

We refer to fll.2l) - fll.3l) and/or (11.41) as the degenerate parametric oscillator [33]. It can be recog- 
nized as a special case of so-called generalized harmonic oscillators that had attracted considerable 
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attention over the years in view of their great importance to many advanced quantum problems 
including coherent states and uncertainty relations, Berry's phase, asymptotic and numerical meth- 
ods, quantization of mechanical systems, Hamiltonian cosmology, charged particle traps and motion 
in uniform magnetic fields, molecular spectroscopy and polyatomic molecules in varying external 
fields. Quadratic Hamiltonians have particular applications in quantum electrodynamics because 
the electromagnetic field can be represented as a set of forced harmonic oscillators. See, for ex- 
ample, [3], [6], [9], [TO], [m, [13], [15], [16], [22], [26], [30], [13], [H], [19] and references therein. 
Nonlinear oscillators play a central role in the theory of Bose-Einstein condensation [8] because the 
dynamics of gases of cooled atoms in a magnetic trap at very low temperatures can be described by 
an effective equation for the condensate wave function known as the Gross-Pitaevskii (or nonlinear 
Schrodinger) equation [IT], [IB], [12], [5i] and [15]. 

A goal of this Letter is to construct Green's function of the Hamiltonian ( ll.4p in the coordinate 
representation in terms of solutions of Ince's equation studied in [25] and [29] (see also references 
therein). Exact time-dependent wave functions and connections with linear and quadratic dynamical 
invariants and SU (1, 1) group are also briefly discussed. 



2. Generalized Harmonic Oscillators 



The degenerate parametric oscillator Hamiltonian f ll.4p belongs to quantum systems described 
by the one- dimensional time- dependent Schrodinger equation 

with variable quadratic Hamiltonians of the form 

d 

H = ait)p^ + hit)x^ + c{t)px + dit)xp, p = -i^^, (2.2) 

ox 

where a{t) , b (t) , c (t) , and d {t) are real- valued functions of time, t, only (see, for example, [1], [S], 
[H], [Z], [3, [n], [n], [13], [23], [2B], [12], [IB] and [ig for a general approach and known elementary 
solutions; a case related to Airy functions is discussed in [21] and our note deals with another special 
case of transcendental solutions). 

The corresponding Green function, or Feynman's propagator, can be found as follows [1], [12]: 

^ = G{x, y, t) = ^ ^i{ao{t)xHMt)-y+io(t)y') ^ (2.3) 



where 



«o(t) = -^^-4T, (2-4) 



4a (t) /io (t) 2a (t) 
X{t) 



/3o (t) = TTT, A (t) = exp / (c (s) - d (s)) ds , (2.5) 

^0 {t) \Jo J 

^°^'^-2/.,(0)/.o(t)^2a(0) ^'-'^ 
and /ig {t) and //^ (t) are two linearly independent solutions of the characteristic equation 

fx" - T (t) /i' + 4(T (t) /X = (2.7) 
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with 



r{t) = — + 2c-2d, a{t) = ab-cd+^l---] (2.8) 



subject to initial conditions 

/io(0) = 0, /i'o(0) = 2a(0)^0, (2.9) 
/ii(0)^0, fi'AO) = 0. (2.10) 

In this note, we present the time-dependent coefficient 7q (t) in terms of two standard solutions /Iq 
and of the characteristic equation ( I2.7p - fl2.8p with initial data f l2.9p - fl2.10p . respectively, instead 
of a previously used formula in terms of only one solution /Xq (see, for example, |6] and |13]). 
Our expression (12. 6 p can be verified by a direct differentiation: 

/^A ' _ /^>o - /^i/^o _ W (/io, /^i) 11) 



2 2 ' 

where the Wronskian can be found with the help of Abel's theorem as follows 

W {fiQ, /^i) = constant (t) a (t) . (2.12) 

Then, as required, 

"^^^ +a{t)Pl{t) = (2.13) 



dt 

and the constant term in (12.60 gives a correct asymptotic as t — )■ 0"^ (cf. Refs. [7j and [12]). The 
Green function f ]2.3p is an eigenfunction of the linear dynamical invariant of Dodonov, Malkin, 
Man'ko and Trifonov [H], [ID], [2S] and (see also Theorem 1 of [15]). 

By the superposition principle the solution of the Cauchy initial value problem can be presented 
in an integral form 

/oo 
G{x,y,t) x{y) dy, \im ij{x,t) = x{x) (2.14) 
t->o+ 



for a suitable initial function x on M (a rigorous proof is given in [42] and uniqueness is analyzed 
in 0; another form of solution is provided by an eigenfunction expansion ^43j; generalized coherent 
states and transition amplitudes are discussed in a forthcoming paper |11]). In the next sections, 
we apply these general results to the case of the degenerate parametric oscillator (11.41) . 

3. Green's Function and Ince Equation 

For the Hamiltonian (11.40 . variable coefficients are 

a = — I 1 + -cos(2a;t) ) , (3.1) 
2m \ OJ J 

mw^ / A , A , , 

c = d = — sin (2ut) (3-3) 
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and general expressions for Green's function fl2.3p -f l2?6|) can be simplified by letting A (t) = 1 and 
c (0) = 0. The characteristic equation is given by 

2Xlo sin (2Lot) , w (w^ _ 3^2^ _ ^ ^^2 ^ ^2ujt) , ^ 

^ + TTTTT^ + ^ = °' (2-^) 

w + A cos (zwrj a; + A cos (zwrj 

which can be identified as a special case of the Ince equation [25], [29] : 

(1 + ao cos 2s) ?/" (s) + 6o sin 2s (s) + (cq + (io cos 2s) y (s) = 0, (3.5) 

when s = ut and parameters are given by 

A A A^ A / A^ \ 
flo = -, bo = 2ao = 2-, Cq = 1 - 3^, do = 1 + — . (3.6) 

00 00 CU^ 00 \ OO"^ J 

Traditionally, a special question which arises in the theory of Ince's equation is the problem of 
the existence of periodic solutions. By Theorem 7.1 on p. 93 of [25], if Ince's equation (13. 5 p has two 
linearly independent solutions of period vr, then the polynomial 

P{0 = 2aoe-M-do/2 (3.7) 

has a zero at one of the points ^ = 0, ±1, ±2, ... . If (13. 5p has two linearly independent solutions of 
period 27i, then 

Q(0 = 2P(e-l/2) (3.8) 
vanishes for one of the values ^ = 0, ±1, ±2, ... . 

In the case of the degenerate parametric oscillator Hamiltonian (13.61) . both of these polynomials 
are strictly positive: 

and 



4A ,.9 A^ 



Thus the corresponding Ince's equation (13. 4p may not have two periodic solutions of periods n/oo 
or 271 /u and our standard solutions /ig and fii that satisfy initial conditions (l2.9l) -( l2TT0l) cannot be 
constructed in terms of simple variants of the Fourier sine and cosine series: 

oo oo 

A^i (^) = ^2n COS {2oont) , /Xq (t) = -B2n sin {2oont) , (3-11) 

71=0 n=l 

or 

oo oo 

/xi(t) = J]A2n+icos(w(2n + l)t), /io W = 1] 52„+i sin (w (2n + 1) t) , (3.12) 

n=0 n=0 

when coefficients and 5^ satisfy certain three-term recurrence relations (more details can be 
found in [25j). Therefore, the degenerate parametric oscillator (11.41) motivates a detailed study of 
non-periodic solutions of Ince's equation. 
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4. The Dynamical Invariant and Time-Dependent Wave Functions 



Exact time-dependent wave functions of generalized harmonic oscillators fl2.1l) - fl2.2p can also 
be obtained as eigenfunctions of the quadratic invariant in terms Hermite polynomials and certain 
solutions of Ermakov-type equation (see Ref. [l3j for a modern introduction and references therein). 
The dynamical invariant of the Hamiltonian f ll.4p can be derived from a general form presented in 
[6], 113] as follows: 

m X sin {2ijjt) — fi' \^ Cq 2 ■ ^ 



\ n 1 + [X/uj) cos {2ujt) J /i^ ox 

where Co > is a constant and function /i (t) satisfies the nonlinear auxiliary equation: 

„ 2Xu sin {2ujt) , u {u^ - 3A^) - X {u^ + X^) cos {2ut) 
^ u + Xcos{2ujt)^ u + Xcos{2ut) ^ 

^(nV {l + {X/u)cos{2ujt)f _ 

A general solution of this nonlinear equation can be found by the following 'law of cosines': 

/i^ it) = Au^ (t) + Bv^ (t) + 2Cu (t) V (t) (4.3) 

in terms of two linearly independent solutions u and v of the homogeneous equation. The constant 
Cq is related to the Wronskian of two linearly independent solutions u and v : 

AB-C' = CojJ^^,, Wiu,v) = uv'-u'v (4.4) 

(more details can be found in [13]; see also references therein regarding this so-called Finney's 
solution). Thus solution of the corresponding initial value problem is given by 

/"'(o) , /i(o) , Co 2 



= U^-^o W + ^-"1 W +^/^ -"(0)^0 (4.5) 

\2a{0) /ii(0) J /i2(0) 

in terms of standard solutions /ig and fii corresponding to initial data fl2.9p -( ETTU]) of the homoge- 
neous equation (14.21) . 

Time-dependent wave functions can be presented in the form 

{x, t) = e-^("+i/2)vW (^^ t) (4.6) 
(see, for example, [13], [19] and references therein), where 

at m jj,'^ 

and orthonormal time-dependent eigenfunctions {x-, t) of the quadratic invariant ( 14. ip are ex- 
pressed in terms of Hermite polynomials |32j : 

{x,t) = Z}„exp zx^— 7^7;, , 4.8 

\ 2n 1 + (A/a;) cos (2a;t) / 



/i / ' " A/7r2"n!/i 
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Then 



^^ = HiJ^, Eij^ = 2^o(n+^]i^^. (4.9) 



If the orthonormal initial wave function is given by 



(x, 0) = e(^^- /')^ (ex) , (4.10) 

where e and 6 are constants, then initial conditions for the auxiliary equation (14. 2 p are given by 

M0) = ^, /.'(0)=2Co^/^A^ + AV (4.n) 
e m \ uj J e 

(the integral Co > can be chosen at reader's convenience). Then integration of Ermakov's equation 
(14. 2p . say with the help of Pinney's solution (14. Sp . determines complete dynamics of the oscillator- 
type wave function through the explicit formula (14. 8p . Further relations of these particular solutions 
of the Schrodinger equation with the Cauchy initial value problem are discussed in [43] (see Theo- 
rem 2). 

5. The SU (1,1) Symmetry 

The Hamiltonian (ll.2p can also be rewritten as follows 

H (t) = JkuJo - n\ (e^*"* J+ + J_) , (5.1) 
where generators of a non-compact SU (1, 1) group are given by 

Then a use can be made of the group properties of the corresponding discrete positive series P;^ 
for further investigation of the degenerate parametric oscillator. This is a 'standard procedure' — 
more details can be found in Refs. [23], [26], [28], [32] and [IQ] and/or elsewhere. 

6. A Special Case 

An integrable special case A = a; = m = ^= lof the degenerate parametric oscillator has been 
recently considered by Meiler, Cordero-Soto and Suslov [28], [7] (this Hamiltonian is a simplest 
time-dependent quadratic integral of motion for the linear harmonic oscillator [6]). Here, the Ince 
equation (13. 4 p simplifies to 

/i" + 2tant/x'- 2/2 = 0. (6.1) 

It has two elementary non-periodic standard solutions: 

/Xq = cost sinh t -|- sin t cosh t, /Xq (0) = 0, /Xq (0) = 2, (6.2) 
/i]^ = cost cosh t -|- sin t sinh t, /i^ (0) = 1, /i'^ (0) = (6.3) 

and Green's function is given in terms of trigonometric and hyperbolic functions as follows 

G{x,y,t) = ^ = (6.4) 

^y 2Txi (cos t sinh t + sin t cosh t) 

(x^ — y"^) sin t sinh t + 2xy — (x^ + y"^) cos t cosh t ' 



X exp 



2i (cos t sinh t -|- sin t cosh t) 
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as a simple consequence of expressions fl2.3p - fl2.6p . It is worth noting that formula (16 ■4p has been 
obtained in [2H] by a totally different approach using the SU (1, l)-symmetry of n-dimensional 
oscillator wave functions and properties of the Meixner-Pollaczek polynomials. More details can 
be found in [28], [7]. 



7. An Extension 



The degenerate parametric amplification with time-dependent amplitude and phase had been 
considered by Raiford [38]. The corresponding Hamiltonian, without damping and neglecting high- 
frequency terms, has the form 



H{t) = -^ {a at + at a) - (^e'(2-<+5(*)) ^2 ^ ^-^(2.t+sit)) ^^iyj ^^ -^-^ 

In this model, the phenomenological coupling parameter A (t) , which describes the strength of the 
interaction between the quantized signal of frequency u and the classical pump of frequency 2u, 
and the pump phase S (t) are in general functions of time as indicated. It includes the special case 
of the pump and signal being off- resonance by a given amount e, i. e., the pump frequency being 
2u + e, by letting 6 (t) = et and A {t) = A, a constant. More details, including Heisenberg's equation 
of motion for the annihilation and creation operators for the signal photons, can be found in the 
original paper [38] . 

In the coordinate representation, the Hamiltonian (17. ip is rewritten as 

H = — (l + ^cos{2ut + 5{t))]p' (7.2) 



2m \ oj 

mco'^ ( ^ \(t) . r / xx\ 9 

+ i 1 - cos {2ut + 5 (t)) j 

+ ^ sin (2ujt + 6 (t)) (px + xp) , P = -S- 
2 t ox 

and Green's function can be found by using the method described in section 2. Further details are 

left to the reader. 



8. A Conclusion 



In this Letter, we have constructed Green's function of the degenerate parametric oscillator (ll.4p 
in terms of standard non-periodic solutions of Ince's equation. To the best of our knowledge, this 
oscillator was introduced for the first time by Takahasi [12] in order to describe the process of 
degenerate parametric amplification in quantum optics (see also [21], [33], [35], [37], [3H], [S], [E] 
and references therein). The corresponding Hamiltonian had also been considered later by Angelow 
and Trifonov [2] , [1] in order to describe the light propagation in a nonlinear anisotropic waveguide. 
Our observation motivates further investigation of properties of the degenerate parametric oscillator 
including a systematic study of corresponding non-periodic solutions of Ince's equation that seems 
to be missing in mathematical literature. Moreover, the Hamiltonian (17. ip . which was considered by 
Raiford [38j in the case of the degenerate parametric amplification with time- dependent amplitude 
and phase, requires a generalization of Ince's equation. 

Acknowledgments. We thank Professor Carlos Castillo-Chavez, Professor Vladimir I. Man'ko, 
Professor Svetlana Roudenko and Professor Dimiter A. Trifonov for support, valuable discussions 



8 



RICARDO CORDERO-SOTO AND SERGEI K. SUSLOV 



and encouragement. The authors are indebted to Professor Victor V. Dodonov for kindly point- 
ing out papers [13], [2H] and |3T] on the degenerate parametric amphfication to our attention — 
his valuable comments have essentially helped to improve the paper. This work is done as a 
part of the summer 2010 program on analysis of Mathematical and Theoretical Biology Institute 
(MTBI) and Mathematical, Computational and Modeling Sciences Center (MCMSC) at Arizona 
State University. The MTBI/SUMS Summer Research Program is supported by The National Sci- 
ence Foundation (DMS-0502349), The National Security Agency (dod-h982300710096). The Sloan 
Foundation and Arizona State University. One of the authors (RCS) is supported by the following 
National Science Foundation programs: Louis Stokes Alliances for Minority Participation (LSAMP): 
NSF Cooperative Agreement No. HRD-0602425 (WAESO LSAMP Phase IV); Alliances for Grad- 
uate Education and the Professoriate (AGEP): NSF Cooperative Agreement No. HRD-0450137 
(MGEOMSA AGEP Phase II). 

References 

[1] A. Angelow, Light propagation in nonlinear waveguide and classical two-dimensional oscillator, Physica A 256 
(1998) #7, 485-498. 

[2] A. Angelow and D. A. Trifonov, Schrddinger covariance states in anisotropic waveguides, Preprint ICTP, 

IC/95/44 Trieste, Italy, 1995, 15 pp. 
[3] M. V. Berry, Classical adiabatic angles and quantum adiabatic phase, J. Phys. A: Math. Gen. 18 (1985) # 1, 

15-27. 

[4] R. Cordero-Soto, R. M. Lopez, E. Suazo and S. K. Suslov, Propagator of a charged particle with a spin in 

uniform magnetic and perpendicular electric fields, Lett. Math. Phys. 84 (2008) #2-3, 159-178. 
[5] R. Cordero-Soto, E. Suazo and S. K. Suslov, Models of damped oscillators in quantum mechanics , Journal of 

Physical Mathematics, 1 (2009), S090603 (16 pages). 
[6] R. Cordero-Soto, E. Suazo and S. K. Suslov, Quantum integrals of motion for variable quadratic Hamiltonians, 

Annals of Physics, article in press (2010), doi:10.1016/j.aop.2010.02.020; see also iarXiv:0912.490 0^^9 [math-ph] 

19 Mar 2010. 

[7] R. Cordero-Soto and S. K. Suslov, Time reversal for modified oscillators , Theoretical and Mathematical Physics 

162 (2010) #3, 286-316; see also arXiv:0808.3149i^9 [math-ph] 8 Mar 2009. 
[8] F. Dalfovo, S. Giorgini, L. P. Pitaevskii and S. Stringari, Theory of Bose-Einstein condensation in trapped gases. 

Rev. Mod. Phys. 71 (1999), 463-512. 
[9] V. V. Dodonov, I. A. Malkin and V. I. Man'ko, Integrals of motion, Green functions, and coherent states of 

dynamical systems. Int. J. Theor. Phys. 14 (1975) # 1, 37-54. 
[10] V. V. Dodonov and V. I. Man'ko, Invariants and correlated states of nonstationary quantum systems, in: In- 
variants and the Evolution of Nonstationary Quantum Systems, Proceedings of Lebedev Physics Institute, vol. 

183, pp. 71-181, Nauka, Moscow, 1987 [in Russian]; English translation published by Nova Science, Commack, 

New York, 1989, pp. 103-261. 
[11] R. P. Feynman, The Principle of Least Action in Quantum Mechanics , Ph. D. thesis, Princeton University, 

1942; reprinted in: "Feynman's Thesis - A New Approach to Quantum Theory", (L. M. Brown, Editor), World 

Scientific Publishers, Singapore, 2005, pp. 1-69. 
[12] R. P. Feynman, Space-time approach to non-relativistic quantum mechanics. Rev. Mod. Phys. 20 (1948) # 2, 

367-387; reprinted in: "Feynman's Thesis - A New Approach to Quantum Theory", (L. M. Brown, Editor), 

World Scientific Publishers, Singapore, 2005, pp. 71-112. 
[13] R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals, McGraw-Hill, New York, 1965. 
[14] S. Fliigge, Practical Quantum Mechanics, Springer- Verlag, Berlin, 1999. 

[15] D. Gomez Vergel and E. J. S. Villasehor, The time- dependent quantum harmonic oscillator revisited: Applications 

to quantum theory, Ann. Phys. 324 (2009), 1360-1385. 
[16] J. H. Hannay, Angle variable holonomy in adiabatic excursion of an integrable Hamiltonian, J. Phys. A: Math. 

Gen 18 (1985) # 2, 221-230. 
[17] Yu. Kagan, E. L. Surkov and G. V. Shlyapnikov, Evolution of Bose- condensed gas under variations of the 

confining potential, Phys. Rev. A 54 (1996) #3, R1753-R1756. 



THE DEGENERATE PARAMETRIC OSCILLATOR 



9 



[18] Yu. Kagan, E. L. Surkov and G. V. Shlyapnikov, Evolution of Bose gas in anisotropic time- dependent traps, 

Phys. Rev. A 55 (1997) #1, R18-R21. 
[19] Yu. S. Kivshar, T. J. Alexander and S. K. Turitsyn, Nonlinear modes of a macroscopic quantum oscillator, 

Phys. Lett. A 278 (2001) #1, 225-230. 
[20] L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Nonrelativistic Theory, Pergamon Press, Oxford, 1977. 
[21] N. Lanfear and S. K. Suslov, The time- dependent Schrddinger equation, Riccati equation and Airy functions, 

[irXiv:0903.3608i^5 [math-ph] 22 Apr 2009. 
[22] P. G. L. Leach, Berry's phase and wave functions for time- dependent Hamiltonian systems, J. Phys. A: Math. 

Gen 23 (1990), 2695-2699. 

[23] C. F. Lo, Coherent-state propagator of the generalized time- dependent parametric oscillator, Europhys. Lett. 24 
(1993) #5, 319-323. 

[24] W. H. Louiseh, Radiation and Noise in Quantum Electronics, McGraw-HiU, Inc., New York, 1964. 
[25] W. Magnus and S. Winkler, Hill's Equation, Dover Publications, New York, 1966. 

[26] I. A. Malkin and V. I. Man'ko, Dynamical Symmetries and Coherent States of Quantum System, Nauka, Moscow, 
1979 [in Russian]. 

[27] I. A. Malkin, V. I. Man'ko and D. A. Trifonov, Linear adiabatic invariants and coherent states, J. Math. Phys. 
14 (1973) #5, 576-582. 

[28] M. Meiler, R. Cordero-Soto, and S. K. Suslov, Solution of the Cauchy problem for a time- dependent Schrddinger 
equation, J. Math. Phys. 49 (2008) #7, 072102: 1-27; see also arXiv: 0711.0559v4 [math-ph] 5 Dec 2007. 

[29] R. Mennicken, On Ince's equation. Archive for Rational Mechanics and Analysis 29 (1968) #2, 144-160. 

[30] S. Menouar, M. Maamache and J. R. Choi, An alternative approach to exact wave functions for time- dependent 
coupled oscillator model of charged particle in variable magnetic field. Annals of Physics, article in press (2010), 
doi:10.1016/j.aop.2010.04.011. 

[31] E. Merzbacher, Quantum Mechanics, third edition, John Wiley & Sons, New York, 1998. 

[32] A. F. Nikiforov, S. K. Suslov, and V. B. Uvarov, Classical Orthogonal Polynomials of a Discrete Variable, 

Springer- Verlag, Berlin, New York, 1991. 
[33] M. Orszag, Quantum Optics, Springer- Verlag, Heidelberg, 1994. 

[34] V. M. Perez-Garcia, P. Torres and G. D. Montesinos, The method of moments for nonlinear Schrddinger equa- 
tions: theory and applications, SIAM J. Appl. Math. 67 (2007) #4, 990-1015. 

[35] J. Pefina, Z. Hradil, and B. Jurco, Quantum Optics and Fundamentals of Physics, Kluwer, Dordrecht, 1994. 

[36] R. R. Puri, Mathematical Methods of Quantum Optics, Springer- Verlag, Berlin, Heidelberg, 2001. 

[37] M. T. Raiford, Statistical dynamics of quantum oscillators and parametric amplification in a single mode, Phys. 
Rev. A 2 (1970) #4, 1541-1558. 

[38] M. T. Raiford, Degenerate parametric amplification with time- dependent pump amplitude and phase, Phys. Rev. 
A 9 (1974) #5, 2060-2069. 

[39] Y. R. Shen, Principles of Quantum Optics, Wiley, New York, 1984. 

[40] Yu. F. Smirnov and K. V. Shitikova, The method of K harmonics and the shell model, Soviet Journal of Particles 

& Nuclei 8 (1977) #4, 344-370. 
[41] D. Stoler, Photon antibunching and possible ways to observe it, Phys. Rev. Lett. 33 (1974) ^^23, 1397-1400. 
[42] E. Suazo and S. K. Suslov, Cauchy problem for Schrddinger equation with variable quadratic Hamiltonians , 

under preparation. 

[43] S. K. Suslov, Dynamical invariants for variable quadratic Hamiltonians , Physica Scripta 81 (2010) #5, 055006 

(11 pp); see also larXiv:1002.0144l ^6 [math-ph] 11 Mar 2010. 
[44] S. K. Suslov, Coherent states and transition amplitudes for generalized harmonic oscillators , under preparation. 
[45] H. Takahasi, Information theory of quantum- mechanical channels, in: Advances in Communication Systems. 

Theory and Applications, v. 1, ed (A. V. Balakrishnan, Editor), Academic Publishers, New York, pp. 227-310. 
[46] C. Trallero-Giner, J. Drake, V. Lopez-Richard, C. Trallero-Herrero and J. L. Birman, Bose-Eistein condensates: 

Analytical methods for the Cross-Pitaevskii equation, Phys. Lett. A 354 (2006), 115-118. 
[47] D. F. Walls and G. J. Milburn, Quantum Optics, Springer- Verlag, Berlin, Heidelberg, 1994. 
[48] K. B. Wolf, On time- dependent quadratic Hamiltonians, SIAM J. Appl. Math. 40 (1981) #3, 419-431. 
[49] K-H. Yeon, K-K. Lee, Ch-I. Um, T. F. George and L. N. Pandey, Exact quantum theory of a time- dependent 

bound Hamiltonian systems, Phys. Rev. A 48 (1993) # 4, 2716-2720. 



10 RICARDO CORDERO-SOTO AND SERGEI K. SUSLOV 

Mathematical, Computational and Modeling Sciences Center, Arizona State University, Tempe, 
AZ 85287-1904, U.S.A. 

E-mail address: ricardojavierSlQgmail.com 

Mathematical, Computational and Modeling Sciences Center, Arizona State University, Tempe, 
AZ 85287-1904, U.S.A. 

E-mail address: sks@asu.edu 

URL: http : //halm . la . asu . edu/ ~ suslov/ index . html 



